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Elemental 2D materials exhibit intriguing heat transport and phononic properties. Here we have 
investigated the lattice thermal conductivity of newly proposed arsenene, the 2D honeycomb struc¬ 
ture of arsenic, using ab initio calculations. Solving the Boltzmann transport equation for phonons, 
we predict a highly anisotropic thermal conductivity, of 30.4 and 7.8 W/mK along the zigzag and 
armchair directions, respectively at room temperature. Our calculations reveal that phonons with 
mean free paths between 20 nm and 1 pm provide the main contribution to the large thermal con¬ 
ductivity in the zig-zag direction, mean free paths of phonons contributing to heat transport in 
the armchair directions range between 20 and 100 nm. The obtained low and anisotropic thermal 
conductivity, and feasibility of synthesis, in addition to other reports on high electron mobility, 
make arsenene a promising material for a variety of applications, including thermal management 
and thermoelectric devices. 

PACS numbers: 


The discovery of graphene as a stable atomically thin 
material has led to extensive investigation of similar 2D 
systems. Its properties such as high electron mobility^, 
and very high thermal conductivity^Ti^ make graphene 
very appealing for applications in electronics, for packag¬ 
ing and thermal management^Ai. The successful iso¬ 
lation of single-layer graphene fostered the search for 
further ultra-thin 2D structures, such as silicene, ger- 
manene, phosphorene, and transition metal dichalco- 
genides, e.g. M 0 S 2 and WS 2 ^i^. These materials are 
now considered for various practical usages due to their 
distinguished properties stemming from their low dimen¬ 
sionality. Thermal transport in 2D materials has recently 
attracted the attention of the scientific community, as 
anomalous heat conduction has been predicted to oc¬ 
cur in systems with reduced dimensionalityi^. Phononic 
properties and thermal conductivity vary significantly 
from one 2D system to another— For example, sil¬ 
icene has a buckled structure and a lower thermal con¬ 
duct ivityiSii^ compared to graphen o^^’^^i^^ . 

2D structures of arsenic and phosphorous have been 
recently investigated^^— . Arsenic and phosphorus are 
in the 5th group of the periodic table and both have dif¬ 
ferent allotropes. Black phosphorus is a layered allotrope 
of phosphorus similar to graphite, and the stability of its 
single layer form, named phosphorene has been probed 
both theoretically and experimentall y . Gray arsenic 
is one of the most stable allotropes of arsenic with a 
buckled layered structure^li^. In addition, arsenic has 
an orthorhombic phase (puckered) similar to black phos- 
phorus2^i2^i2£, and its monolayer is called arsenene (see 
Fig.[T]). Experimental observations have shown that gray 


arsenic undergoes a structural phase transition to the 
orthorhombic precursor of arsenene at temperatures of 
about T = 370 K^. As a monolayer arsenene has a direct 
band gap as opposed to the multilayer allotrope, which 
exhibits an indirect band gap of the order of 1 eV— . 
According to our calculations, monolayer of arsenene is 
stable also near zero temperature, in agreement with pre¬ 
vious reports^^i^i^. Here we investigate heat conduc¬ 
tion in arsenene and we elucidate the anisotropy of its 
thermal conductivity. Since in semiconductors phonons 
are the predominant heat carriers, we use first-principles 
anharmonic lattice dynamics calculations and the Boltz¬ 
mann transport equation to compute phonon dispersion 
relations and thermal conductivity. 

We employ density functional theory as imple¬ 
mented in the VASP code^ with projector-augmented- 
wave (PAW)^ pseudopotentials and the Perdew-Burke- 
Ernzerhof (PBE) generalized gradient functional — 
exchange-correlation functional. Integration over Bril- 
louin zone is performed using 15 x 15 x 1 F-center 
Monkhorst-Pack mesh of k-points^i and plane waves up 
to an energy cutoff of 360 eV are used as basis set to 
represent the Kohn-Sham wave functions. The simula¬ 
tion cell is built with a vacuum layer of 10 A, which is 
sufficient to avoid interactions between periodic images. 
The structural relaxations were performed by using the 
conjugate-gradient (CG) algorithm with a convergence 
criterion of 10“® eV /A for the maximum residual force 
component per atom. 

As shown in Fig. [I] arsenene has a puckered structure 
arranged in a rectangular unit cell with lattice vectors 
a = 3.686 A and b = 4.769 A. The bond length between 
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FIG. 1: Top and side view of puckered arsenene. 



FIG. 2: (a) Phonon dispersion curves along the high- 

symmetry directions of the first Brillouin zone for puckered 
arsenene. 


TABLE I: The lattice parameters of arsenene. a and b are 
the lattice vectors, A is the puckering distance, ci is the bond 
length between atoms that are in the same plane, C 2 is the 
bond length between atoms that are in different planes. 

Mono-layer of arsenene 
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atoms in the same plane (upper panel or lower) is 2.510 
A, while the bonding length between atoms in different 
planes is 2.496 A. The distance between the two planes 
is denoted by A and is 2.402 A. Our structural param¬ 
eters are compared with previous studies in Table H We 
have utilized the PHONOPY package^ interfaced with 
VASP to compute inter-atomic force constants (IFCs). 
After testing convergence of the phonon frequencies as a 
function of the size of the supercell, a 8 x 8 x 1 super-cell 
has been adopted, with a 3 x 3 x 1 P-center Monkhorst- 
Pack mesh of k-points, to build the dynamical matrix T>. 
The frequencies of the vibrational modes are obtained by 
diagonalizing T), and the phonon dispersion curves are 
represented in Fig. [21 The absence of phonon branches 
with negative phonon frequencies indicates the stability 
of the system at zero temperature. 

The phonon bands exhibit two acoustic modes, a lon¬ 
gitudinal one (LA) and an in-plane transverse one (TA), 
with linear dispersion (oj oc q) for q — t 0, and an 
out-of-plane flexure mode (ZA) with quadratic disper¬ 
sion in the long wavelength limit, which is a general 
feature of 2D membranes^^ and have been characterized 
in graphenei^i^l, silicen o^^i^^ , hBN^, MoS 2 ^ and ultra 
thin silicon membranes^i^. Flexural modes have a ma¬ 
jor role in contributing to the thermal conductivity of 
graphene, both as carriersi^ and as scatterers^. In ad¬ 
dition, we observe a gap of about 2 THz between the 
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FIG. 3: The thermal conductivity of arsenene as a function 
of temperature along the zigzag and armchair direction. 


low-frequency and high-frequency optical branches. 

We have computed the thermal conductivity of ar¬ 
senene by solving the phonon Boltzmann transport 
equation (BTE) with an iterative self-consistent algo- 
rith m^'^d^ , using the ShengBTE package^. In our calcu¬ 
lations we have considered third-order interatomic force 
constants up to the forth shell of neighbors in a 5 x 5 x 1 
super-cell with a 3 x 3 x 1, F-centered Monkhorst-Pack 
mesh of k-points^i^. Intrinsic three phonon scatter¬ 
ing processes have been considered in the calculation of 
phonon relaxation times: 


.^0 - /V ( X! ^AA'A" + 9 X! 


Errv 




A'A'^ 


( 1 ) 

where and represent the scattering rates 

due to absorbing and emitting three-phonon processes. 
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TABLE II: Sound velocity (group velocity) of acoustic modes 
of arsenene close to the zone center (L-point). (*) Due to 
quadratic form of ZA mode near T point, the magnitude has 
been averaged over a small reigon interval. 


Sound velocity km/s 



Arsenene 

Phosphorene 


ZA (*) 

TA 

LA 

LA (ref.'^) 

Armchair 

0.72 

2.41 

2.44 

4.17 

Zigzag 

0.34 

2.41 

4.77 

7.73 


respectively. The extrinsic scattering term is zero 
in our calculations, since arsenic does not present iso¬ 
topic disorder, and neither defect or boundary scatter¬ 
ing processes were considered. Knowing the distribution 
function, the phononic thermal conductivity tensor (k) 
is expressed as: 

> ( 2 ) 

where / is the phonon distribution function, uj\ is the 
phonon angular frequency, is the velocity component 
along the a-direction, and tx is the phonon relaxation 
time. 

The two diagonal components of the thermal conduc¬ 
tivity at temperatures between 20 and 500 K are re¬ 
ported in Fig. |3l The thermal conductivity along the 
zigzag and armchair direction at T = 300 K is 30.4 and 
7.8 Wm“^K“^, respectively. The nearly 3-fold ratio be¬ 
tween the two components indicates a high anisotropy. 
The thermal conductivity of aresenene in comparison to 
graphene and hexagonal boron nitride is very low: for 
example both experimental and theoretical studies in¬ 
dicate the intrinsic thermal conductivity of graphene is 
larger than 3000 W/mK — Due to the larger atomic 
weight of As, k of arsenene is also lower than that of 
phosphorene, which is also strongly anisotropic, theoret¬ 
ically predicted as 110 and 36 Wm“^K“^ in the zigzag 
and armchair directions, respectively^^. We note indeed 
that the group velocity of the acoustic modes of arsenene 
along both zigzag and armchair directions near the zone 
center are about twice as smaller as those of phosphorene 
(see Table HB. Since according to Equation [5] group ve¬ 
locities contribute to n quadratically, we can argue that 
the smaller vx is responsible for the difference between 
the thermal conductivities of arsenene and phosphorene. 

The contribution of each phonon mode to the total 
thermal conductivity as a function of temperature is 
shown in Fig. |4l The ZA mode provides the largest 
contribution to thermal conductivity along both direc¬ 
tions at low temperatures. However, this contribution 
decreases as the temperature increases and reaches 7.8% 
(armchair direction) for temperatures above 300 K. This 
is in contrast to the case of graphene where the ZA modes 
provide the main contribution to heat conductioni^. Due 
to the special selection rules of phonon scattering in 
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FIG. 4: Normalized contribution of each phonon branch to 
the total thermal conductivity, as a function of temperature 
along the (a) zigzag and (b) armchair directions. 


graphene, the ZA mode has a relatively large relaxation 
time compared to other phonon modes^^, while due to 
the puckered structure of arsenene, similar to silicene and 
phosphoren o^^d^ , the hexagonal symmetry of graphene 
is broken, resulting in larger scattering rates and shorter 
relaxation times of the ZA modes at high temperatures. 
The contribution of the ZA modes remains however sub¬ 
stantial and roughly equivalent to that of the TA modes 
in the zigzag direction, while it drops below 10% in the 
armchair direction. 

At temperatures higher than 100 K, the contribution 
of LA modes to the thermal conductivity is the largest in 
both directions. The group velocity of LA mode near the 
zone center is the largest along the zigzag direction, while 
it is comparable to the TA mode along the armchair di¬ 
rection (see Table |ll| . Furthermore, the relaxation times 
of the LA modes at low frequencies, which play a signif¬ 
icant role in the process of heat conduction, are larger 
than that of other acoustic modes (see Fig.[S]). All these 
effects combined together result in a larger contribution 
of LA modes to the thermal conductivity at high tem¬ 
peratures. 

Finally we computed the cumulative thermal conduc¬ 
tivity as a function of phonon mean free path, defined 
as 

ki{A ^ Aniax) = "" An > (3) 

f I A 

^^vn.ax 

where Amax is the maximum mean free path and Aq is 
a fitting parameter that can be interpreted as the mean 
free path of heat-carrying phonons of infinite system size. 
By fitting Eq.(j3l) to the calculated values one can eval¬ 
uate Aq, see Eig. [S] This quantity shows which mean 
free paths contribute the most to the thermal conduc¬ 
tivity. Phonon modes with mean free paths larger than 
59.5 nm along the armchair direction and 91.8 nm along 
the zigzag direction have 50% contribution to the to¬ 
tal thermal conductivity. Hence we can argue that in a 
polycrystalline structure with a grain size of the order of 
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FIG. 5: Total relaxation time of phonon modes as a function 
of frequency v along the (a) zigzag and (b) armchair direction. 



Maximum MFP (nm) 


50 nm the thermal conductivity can be significantly fur¬ 
ther reduced. However, given the broad range of mean 
free paths that contribute, the most efficient approach to 
abate k would be hierarchical nanostructuring^. 


In summary we have theoretically investigated the 
phonon properties and the thermal conductivity of ar- 
senene by ab initio anharmonic lattice dynamics. Our 
calculations evidence the lowest thermal conductivity 
among the 2D elemental materials found so far. The 
larger mass of arsenic atoms compared to phosphorus 
leads to a 3 times smaller thermal conductivity of this 
material compared to phosphorene. Similar to phos- 
phorene, the thermal conductivity of arsene is strongly 
anisotropic. Such anisotropy could be utilized in ther¬ 
mal engineering application, to establish directional heat 
transport at the nanoscale. Unlike graphene, the longitu¬ 
dinal acoustic phonon modes provide the main contribu¬ 
tion to the thermal conductivity of arsenene at tempera¬ 
tures above 100 K, which is mainly due to the puckered 
structure of this material. From this observation, we ex¬ 
pect that K of arsenene would be less affected by the inter¬ 
action with a substrate than in those materials for which 
the main contribution comes from flexural modes. In ad¬ 
dition, given the high electron mobility and the reduced 
dimensionality, which should enhance the Seebeck coeffi¬ 
cient, its relatively low and strongly directional thermal 
conductivity makes arsenene a good candidate from ther¬ 
moelectric energy conversion. 


FIG. 6: Cumulative thermal conductivity of arsenene as a 

function of the phonon mean free path at T = 300 K along The work of S.M.V.A. was supported in part by the 
the zigzag and armchair directions. Research Council of the University of Tehran. 
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